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Abstract. We have assembled a large data-set of 613 
galaxy redshifts in the Coma cluster, the largest presently 
available for a cluster of galaxies. We have defined a sam- 
^ pie of cluster members complete to b26 5 = 20.0, using a 
O membership criterion based on the galaxy velocity, when 
available, or on the galaxy magnitude and colour, oth- 
QQ erwise. Such a data set allows us to define nearly com- 
^— H plete samples within a region of 1 h~^ Mpc radius, with a 
sufficient number of galaxies per sample to make statisti- 
cal analyses possible. Using this sample and the ROSAT 
T-H PSPC X-ray image of the cluster, we have re-analyzed 
the structure and kinematics of Coma, by applying the 
wavelet and adaptive kernel techniques. 

A striking coincidence of features is found in the dis- 
tributions of galaxies and hot intracluster gas. The two 
^ central dominant galaxies, NGC 4874 and NGC 4889, are 
rS |Surrounded by two galaxy groups, mostly populated with 
^ galaxies brighter than b26 5 = 17 and well separated in 
velocity space. On the contrary, the fainter galaxies tend 
to form a single smooth structure with a central peak co- 
inciding in position with a secondary peak detected in X- 
rays, and located between the two dominant galaxies; we 
suggest to identify this structure with the main body of 
the Coma cluster. A continuous velocity gradient is found 
in the central distribution of these faint galaxies, a prob- 
able signature of tidal interactions rather than rotation. 
There is evidence for a bound population of bright galaxies 
around other brightest cluster members. 

Altogether, the Coma cluster structure seems to be 
better traced by the faint galaxy population, the bright 
galaxies being located in subclusters. We discuss this evi- 
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dence in terms of an ongoing accretion of groups onto the 
cluster. 
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1. Introduction 

The Coma cluster of galaxies, number 1656 in the Abell 
(1958) catalogue, has been one of the most studied galaxy 
clusters since the first determination of its mass by Zwicky 
(1933). Its high richness (richness class 2) and smooth ap- 
pearance made of this cluster the prototype of well relaxed 
clusters, in contrast to the most nearby cluster, Virgo, 
considered as the prototype of irregular clusters. Support 
to this view came from the discovery by Rood (1969) of 
luminosity segregation among the galaxies in Coma, inter- 
preted as an indication of partial dynamical relaxation. 

However, as early as in the 70's some facts started 
to challenge the picture of Coma as a well relaxed dy- 
namical system: Rood & Sastry (1971) classified Coma as 
a "binary" cluster, because of the presence of a couple 
of giant galaxies in its core (NGC 4889 and NGC 4874), 
and Rood (1974) hinted to a possible subclustering around 
these two galaxies. Valtonen & Byrd (1979) went even far- 
ther, by proposing an unconventional dynamical model in 
which most of the cluster mass was linked to these two gi- 
ant galaxies. During the 80's the claims for subclustering 
in Coma became more and more frequent: Baier (1984) 
claimed significant substructure in the South- West region 
of Coma, around the giant galaxy NGC 4839, and hinted 
at a possible subclustering in the cluster core; Gerbal et 
al. (1986) and Mellier et al. (1988) suggested that subclus- 
tering occurs around many bright galaxies and is not only 
limited to NGC 4889, 4874 and 4839; Fitchett k Webster 
(1987) used the Lee-statistics to provide the first signifi- 
cant detection of subclustering in the cluster core. 
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Nevertheless, in the same period, Dressier & Sheet- 
man (1988) did not find a significant signal for subcluster- 
ing in Coma, and most dynamical analyses still assumed 
a smooth cluster structure (Kent & Gunn 1982; The & 
White 1986; Merritt 1987; Hughes 1989). It was in fact 
not until the 90 's that the complicated structure of the 
Coma cluster became widely accepted, thanks to new tech- 
niques applied to the galaxy distribution (Escalera et al. 
1992) and to the detection of these structures in the X- 
ray band with the Spacelab 2 X-ray telescope (Watt et 
al. 1992) and particularly with _ROS'^T (Briel et al. 1992; 
White et al. 1993; Vikhlinin et al. 1994). 

Subclustering implies an incomplete virialization of the 
cluster, with ongoing mergers of groups onto the main 
body of the cluster; the dynamical youth of the Coma clus- 
ter is also indicated by the kinematical difference between 
the star-forming galaxies and the old galaxy population 
(Sodre et al. 1989; Zabludoff & Franx 1993; Gavazzi et al. 
1995; Colless & Dunn 1995, hereafter CD95); the existence 
of a population of E-|-A spectrum galaxies likely to have 
suffered a burst of star formation within the last 2 Gyr 
(Caldwell et al. 1993, 1995) is also suggestive of ongoing 
accretion of galaxies onto the cluster. Finally, unlike most 
clusters. Coma lacks a cooling fiow (Edge et al. 1992) and 
has a radio halo (Hanisch 1980); these features have been 
interpreted in terms of an ongoing collision of two equally 
sized subclusters (McGlynn & Fabian 1984; Edge et al. 
1992; Tribble 1993; Fabian et al. 1994). 

Coma can now be considered as the prototype of rich 
clusters endowed with subclusters, and thus not fully re- 
laxed. However, the very presence of substructures log- 
ically implies the existence of a mam body, on which 
these features are superimposed; it has become custom- 
ary to identify this main body with the group surround- 
ing NGC 4874, mainly because this galaxy lies on the main 
peak of the diffuse cluster X-ray emission (see, e.g.. White 
et al. 1993); in particular, this is the point of view of CD95. 
However, the velocity of NGC 4874 is offset from the av- 
erage velocity of the cluster, and although not uncommon 
among brightest cluster members, this is an unexpected 
result if the galaxy is at rest at the bottom of the poten- 
tial well of the cluster. Therefore it is also possible that 
NGC 4874 rather lies at the bottom of the gravitational 
well of lis own group, similarly to what is observed for 
other central dominant galaxies (Bird 1994). 

The large number of galaxy velocities now available 
in the Coma cluster field (613, including our new redshift 
data-set, Biviano et al. 1995a), and the high resolution of 
the ROSAT PSPC X-ray image allows us to examine the 
structure and kinematics of the central part of the Coma 
cluster (within 1 h~^ Mpc, where h is the Hubble constant 
in units of 50 km s~^ Mpc~^), with the aim of understand- 
ing if a mam cluster body actually exists, independently 
from the two central subclusters. In this paper we also 
address the issue of subclustering and ongoing infall onto 
the cluster in a larger cluster region. Some preliminary re- 



sults from our work have been presented in Biviano et al. 
(1994), Gerbal et al. (1994), and Durret et al. (1995). 

Shortly before this paper was submitted, we received 
a preprint from CD95 which also carried out a structural 
analysis of the Coma cluster, with rather different tech- 
niques; while we agree on some of their results, we do not 
come to the same conclusions on what concerns the struc- 
ture of the central cluster region. 

In § 2 we describe the optical and X-ray data used for 
our analysis; we analyse the structure of the central region 
of the Coma cluster, as described by the galaxy and X- 
ray distributions, in § 3, and the cluster kinematics in § 4; 
in § 5 we extend our analysis to the external regions; the 
relevant discussion is provided in § 6; finally we summarize 
our results and give our conclusions in § 7. 

2. The data 

2.1. Optical data 

We have taken positions, blue magnitudes, and blue— red 
colors of galaxies in the Coma region from the catalogue of 
Godwin, Metcalfe k Peach (1983, hereafter GMP). This 
catalogue is nominally complete up to an isophotal mag- 
nitude \)2Gji = 20.0 (we note these isophotal magnitudes 
h hereafter), in a 2.63 square degree area centered on 
ai95o = 12''57™.3, ,^1950 = +28n4'.4. 

We have assembled a total number of 613 galaxies with 
redshifts in the Coma region defined by GMP: redshifts 
for 143 galaxies have been obtained by ourselves through 
observations with the CFHT in May 1993 (Biviano et al. 
1995a); another 167 redshifts have been kindly provided 
by M. Colless in advance of publication (CD95); finally, 
another 303 redshifts have been taken from the literature 
(Bothun et al. 1985; Caldwell et al. 1993; Davies et al. 
1987; da Costa et al. 1984; de Vaucouleurs et al. 1976; 
Gavazzi 1987; Giovanelli & Haynes 1985; Hoffman et al. 
1989; Huchra et al. 1983; Huchra et al. 1990; Huchtmeier 
& Richter 1989; Karachentsev & Kopylov 1990; Kent & 
Gunn 1982; Leech et al. 1988; Lucey et al. 1991; Palumbo 
et al. 1983; Sullivan et al. 1981; van Haarlem et al. 1993; 
Wegner et al. 1990; White et al. 1982; Zabludofl' et al. 
1993). For the compilation of redshifts from the literature, 
we have made use of the ZCAT catalogue (Huchra et al. 
1995). When several redshift measurements were available 
for the same galaxy, we have chosen the most accurate one. 

In Table 1 is indicated the completeness in velocity as a 
function of the limiting magnitude in the GMP region and 
in the inner region (circle of 1500" radius), with respect 
to all galaxies in the field (cols. 2 and 3, respectively). 

For the purpose of our analysis, we have extracted sev- 
eral samples from a combination of the GMP photometric 
catalogue and our redshift catalogue. In order to reach a 
high value of completeness in velocities, and still have a 
sufficiently large number of galaxies in the sample, most of 
our analyses have been performed on a circular region of 
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Table 1. Completeness in velocity as a function of limiting 
magnitude. and N represent the number of galaxies with 
velocities and the total number of galaxies respectively. 



Magnitude range N,j/N 

CMP region 1500" region 



b < 17.0 


282/346 


122/123 


b < 17.5 


325/447 


148/150 


b < 18.0 


387/607 


185/190 


b < 18.5 


453/820 


224/248 


b < 19.0 


508/1128 


259/313 


b < 19.5 


554/1609 


291/403 


b < 20.0 


605/2510 


316/544 


17.0 < b < 19.0 


226/782 


137/190 


17.0 < b < 20.0 


323/2164 


194/421 



radius 1500" (at the distance of the Coma cluster, 1500" 
correspond to 1 h~^ Mpc). We choose a circular region to 
avoid introducing any preferential direction. 

From these samples, we have extracted the cluster 
members. The membership assignment is based on the 
galaxy velocity when available (galaxies with velocities in 
the range 3000 — 10000 kms~^, plus one galaxy with a 
velocity of 2998 km s"-'^, are assumed to belong to Coma), 
or on the location of the galaxy in the colour-magnitude 
band defined by Mazure et al. (1988); for more details on 
this membership assignement see Biviano et al. (1995b). 

We have defined and labeled the samples and subsam- 
ples that we have used as follows: 

— sample C'au contains the galaxies members of the clus- 
ter in the CMP region, up to the completeness magni- 
tude h = 20.0; 

— sample C contains members of the cluster within the 
circle of 1500" radius centered on the CMP center, up 
to & = 20.0; 

— subsample C'iright contains C members with h < 17; 

— subsample C/aj„( contains C members with 17 < & < 
20.0. 

— subsample C/aj„(_i9 contains Cmembers with 17 < 
&<19.0. 

Finally, the subscript "t;" indicates the subsamples of 
galaxies with available velocities. 

The numbers of galaxies with and without velocities 
in these samples are listed in Table 3 (see next section). 

2.2. X-ray data 

X-ray images of Coma obtained with the ROSAT PSPC 
were kindly made available to us by S. White before their 
release in the Garching public archive. We combined two 
of these images, taken in the "hard" energy band (0.4- 
2.4 keV), totalling an exposure time of about 40000 s. The 
image sizes were 256x256 pixels^, with a pixel size of 15" 
(i.e. 10 h~^ kpc at the distance of Coma). 



3. The central region of Coma: structure 

3.1. Methods of analysts 

In order to visualize the distribution of galaxies on all 
scales at once, we have used the method of adaptive ker- 
nels (e.g. Silverman 1986). The adaptive kernel technique 
is widely used in the astrophysical literature (e.g. Pisani 
1993, 1995; Merritt & Gebhardt 1994). Here we have also 
used this technique to draw maps of various moments of 
the galaxy velocity distribution (see Figs. 2 and 12). 

The adaptive kernel technique requires an initial ker- 
nel size; this size is then modified according to the local 
density of points. Unless otherwise stated, we used an ini- 
tial kernel size equal to 200" (i.e. 0.13 h~^ Mpc at the 
distance of the cluster), since this is close to the intrin- 
sic scale of many structures in the Coma cluster (Mellier 
et al. 1988; Escalera et al. 1992), and also happens to be 
nearly equal to the "optimal" kernel size (Silverman 1986) 
for the sample C , which is the most used sample in our 
analysis. 

The significance of the detected structures is estimated 
through a bootstrap resampling analysis. See Appendix A 
for further details. 

For the X-ray data, we have used the wavelet anal- 
ysis, which easily removes noise and allows by its multi- 
scale approach to detect (at a > 3(7 level) both the point- 
sources and the extended sources. A thorough description 
of the method we have used is given in Slezak et al. (1994; 
see also Grebenev et al. 1995), and in Appendix B. The 
wavelet method has already been applied to the investiga- 
tion of both the galaxy and gas distributions in the Coma 
cluster (Escalera et al. 1992; Gerbal et al. 1994; Vikhlinin 
et al. 1994). 

3.2. The galaxy dtstnhutton 

We display in Fig. 1 (upper panel) an adaptive kernel den- 
sity map of all the cluster members in the circle of 1500" 
radius (i.e. sample C ); this and the following adaptive- 
kernel maps, have 128 x 128 pixels. The significance level of 
the features visible in this map has been evaluated through 
a bootstrap resampling technique; the density-map at —3(7 
confidence level (see Appendix A) is also plotted (lower 
panel): most features of the original map are still visible 
in this map. 

A complex structure is observed, with a main peak cor- 
responding to the central dominant galaxy NGC 4874, an 
extension to the East corresponding to the brightest clus- 
ter member, NGC 4889, and a less significant extension to 
the South-East in the direction of another bright galaxy, 
NGC 4911 (these structures have already been pointed out 
by Mellier et al. 1988 and Escalera et al. 1992). This map 
therefore naturally suggests the presence of sub clustering. 

A classical test for estimating the significance level of 
the presence of subclustering is that of Dressier b, Sheet- 
man (1988). In this test one computes the local mean ve- 
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Fig. 1. Upper panel: adaptive kernel density map of all the 
galaxies belonging to the Coma cluster within a circle of ra- 
dius 1500" (sample C). In this figure and in the following ones, 
coordinates are in arcseconds with respect to the cluster cen- 
tre defined in CMP; North is to the top and East to the left; 
isodensity contours are spaced by 1 x 10~^ galaxies/arcsec^ . 
The galaxies NGC 4889, NGC 4874 and NGC 4911 at respec- 
tive positions (-304,22), (124,-41) and (-939,-649) are indicated 
by x's. The position of a secondary maximum detected in the 
X-ray map is indicated by an asterisk at position (-201,-54). 
Lower panel: adaptive kernel density map at —3(7 confidence 
level (see Appendix A). 



locities and velocity dispersions of all possible groups of 11 
neighboring galaxies, and the deviations, of these local 
values from the global cluster mean velocity and velocity 
dispersion (see also Appendix A). Large 6 values indicate 
a high probability for the cluster to contain substructures. 
Bird (1994) suggests that using VWy instead of 11 galax- 
ies in the evaluation of the 6 parameters (where Ny is the 
total number of galaxies with available velocities), max- 
imizes the sensitivity of the test to significant structures 
while reducing sensitivity to noise fiuctuations. 

We have run this test (modified according to the sug- 
gestion of Bird 1994) on our sample Cy , which contains 
255 galaxies with available velocities, and found a prob- 
ability of 0.921 that substructures are present (based 
on 1000 Montecarlo simulations). In order to visualize 
these substructures, we have plotted an adaptive kernel 6- 
map; notice that this is quite comparable to the weighted 
wavelet method developed by Escalera et al. (1994). We 
caution the reader that it is not our intention to improve 
the statistical test for detection of substructure developed 
by Dressier & Shectman (1988); the only purpose of our 
technique is to allow a 2D visualisation of the values of the 
parameter 6 used to identify susbtructures. This is done 
following the procedure described in Appendix A. lso-6 
contours for the sample Cy are shown in Fig. 2. 






Fig. 2. Dressier & Schectman 5-map for sample C,j . The av- 
erage 8 in this map is 40; the first contour starts at a value of 
8 = AO and following contours are spaced by an interval of 20 
in 8. Note that contours beyond a distance of 1500" from the 
center are meaningless, because the C,j sample only contains 
galaxies within this radius. Symbols are the same as in Fig. 1. 



One may see on this figure that there are two main 
peaks in the (5-distribution, corresponding to two regions 
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where galaxy velocities strongly differ from the average 
velocity field; note that contours beyond 1500" from the 
center are meaningless, because the sample only con- 
tains galaxies within this radius. We have selected the 11 
galaxies located within the contours delimiting the two 
main peaks; these galaxies are listed in Table 2. The 6 
galaxies in the first peak have an average velocity of 4518 
± 324 km s~^, and the 5 galaxies in the second peak have 
an average velocity of 8527 ± 265 kms""'^. Notice that 
these two values are close to the boundaries of the range 
of velocities taken for galaxies belonging to Coma. The 
velocity dispersions of these two groups are more in the 
range of those of poor clusters than of groups (see e.g., 
Mazure et al. 1995; Ramella et al. 1995); however, it is 
possible that the first group is in fact made of 3 galaxy 
pairs, and that the galaxy with a velocity of 9811 km s~^ 
is not a member of the second group, but an isolated back- 
ground galaxy. By excluding these 11 galaxies, the test of 
Dressier & Shectman (1988) gives a probability of 0.529 
for the presence of substructures on the remaining 244 
galaxies, and the (5-map does not show any other signif- 
icant feature. So, there are good reasons to believe that 
these groups are in the fore- and background of the Coma 
cluster. In the following, we have excluded these 11 galax- 
ies from our samples. The number of galaxies remaining 
in the subsamples defined in § 2.1 after the elimination of 
these 11 galaxies, are listed in Table 3. 



Table 2. Eleven galaxies belonging to the foreground and 



background 


groups. 








Group 


GMP 


Velocity 


h 


(x,y) 




number 


km 




arcsec 


1 


3262 


3682±100 


16.77 


68,-417 




3275 


3834± 36 


18.66 


76,-419 




3383 


4585± 75 


18.50 


178,-394 




3400 


4692± 16 


15.32 


191,-310 




3473 


5061± 92 


18.92 


249,-409 




3522 


5126± 31 


16.39 


300,-266 


2 


2923 


8664± 51 


17.65 


-302,-711 




3012 


8064± 48 


17.49 


-215,-863 




3071 


8865± 62 


17.17 


-145,-809 




3092 


8434± 68 


17.55 


-126,-629 




3176 


9811± 98 


18.29 


-12,-808 



In Fig. 3 we have displayed the adaptive kernel density 
map and the density map at —3(7 confidence level, for the 
sample C , after the elimination of these 11 galaxies; the 
map appears slightly smoother than before (compare with 
Fig. 1), but the most relevant features are unchanged. In 
particular, the two peaks around the two central dom- 
inant galaxies are still visible. It is not surprising that 
these peaks were not detected by the Dressier & Shect- 
man (1988) test, since the velocity distribution of galaxies 
within these groups is not very different from the overall 
cluster velocity distribution. 



Table 3. Number of cluster members in the various subsam- 
ples, after the elimination of 11 galaxies in two fore- and back- 
ground groups. and N represent the number of galaxies 
with velocities and the total number of galaxies respectively. 



Sample 




Call 


469/1372 


c 


244/372 


Cliright 


118/119 


C f aint 


126/253 


Cfaint,lQ 


108/141 



The structure of Coma as visible from Fig. 3 (upper 
panel) remains complex, despite the elimination of the 
groups of galaxies in the foreground and background, and 
raises the natural question of whether there really is a 
Coma cluster, i.e. an underlying main structure. As we dis- 
cuss in § 6.1, there are reasons to believe that faint galaxies 
may be more faithful tracers of the main cluster body po- 
tential. For this reason, we have split our sample C into the 
two sub-samples C'iright and C faint previously described. 
We have chosen the magnitude 6 = 17 to separate these 
two samples because this magnitude corresponds to the 
dip in the Coma cluster luminosity function (Biviano et 
al. 1995b), and thus naturally separates bright from faint 
galaxies in Coma. Moreover, the number of galaxies with 
available redshifts in the two subsamples is approximately 
the same. 

It turns out that the structures seen in the distribu- 
tion of bright galaxies are quite different from those seen 
in the distribution of faint galaxies, as can be seen by plot- 
ting the adaptive kernel maps of the C'iright (Fig. 4) and 
C faint (Fig. 5) samples respectively. 

As seen in these two figures, the relative importance 
of the group surrounding NGC 4874 is strongly empha- 
sized in the bright galaxy distribution, while the two peaks 
around the two giant central galaxies are no longer present 
in the faint galaxy distribution. The appearance of the 
map shown in Fig. 5 is more regular, with a single den- 
sity maximum located between the position of the two 
central galaxies. This secondary peak is located at GMP 
coordinates (-201,-59), and therefore at distances from 
NGC 4874 and NGC 4889 of 323" and 133" respectively. 
While there still is an overdensity around NGC 4874, it is 
less prominent than before; this overdensity continuously 
decreases (and the secondary peak overdensity increases) 
when the fraction of faint galaxies is increased (i.e. going 
from sample Chright to Cfaint,i9 and to Cfaint )• 

By looking at the density map at —3(7 confidence level 
(lower panel of Fig. 5) one can see that the maximum 
in the faint galaxy distribution has not a very well de- 
termined position and tends to shift towards NGC 4874. 
We have estimated the uncertainty on the position of 
this peak by running 1000 bootstrap resamplings on the 
Cfaint sample, and estimating the position of the highest 
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Fig. 3. Upper panel: adaptive kernel map of all the galaxies 
belonging to the Coma cluster within a circle of radius 1500" 
(sample C), after the elimination of the 11 galaxies belonging 
to the background and foreground groups. The contour spac- 
ing and the symbols are the same as in Fig. 1. Lower panel: 
adaptive kernel density map at —3(7 confidence level (see Ap- 
pendix A). 



peak in each bootstrap sample; the dispersion in the po- 
sitions of the highest peaks is 315". This large dispersion 
is probably due to the presence of two density peaks, the 
one visible in Fig. 5 (upper panel) and the very strong 
peak around NGC 4874 outstanding in the bright galaxy 
distribution and not totally erased in the faint galaxy dis- 
tribution. The striking positional coincidence between the 
density peak of the faint galaxy distribution, and a sec- 



Fig. 4. Upper panel: adaptive kernel map of the galaxies be- 
longing to the Coma cluster with magnitudes b < 17.0 (sample 
Cbr ight )■ The contour spacing and the symbols are the same 
as in Fig. 1. Lower panel: adaptive kernel density map at —3(7 
confidence level (see Appendix A). 



ondary maximum detected in the X-rays (see § 3.3), not 
associated with any of the two central dominant galaxies, 
strongly argues in favour of the reality of this peak. 

Magnitudes therefore appear as an efficient tool of dis- 
crimination: the cluster appearance changes when select- 
ing samples of luminous or faint galaxies. Bright cluster 
galaxies are mainly located around the brightest mem- 
bers, while faint ones delineate a smoother single peaked 
structure. 
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Fig. 5. Upper panel: adaptive kernel map of the galaxies be- 
longing to the Coma cluster with magnitudes b > 17.0 (sample 
C faint)- The contour spacing and the symbols are the same 
as in Fig. 1. Lower panel: adaptive kernel density map at —3(7 
confidence level (see Appendix A). 



3.3. The X-ray dtstnhutton 

We have performed a wavelet analysis on the whole X-ray 
image of Coma, but for the purpose of comparison with 
the optical sample C , we concentrate here only on the 
central part of this image (within a radius of ~ 1500"). 
The reconstructed X-ray image (i.e. that obtained after 
subtracting the noise and artefacts, following a technique 
described in Appendix B) is displayed in Fig. 6. The four 
wavelet images, obtained from the noise-subtracted image 



through the wavelet transform, corresponding to wavelet 
scales of 15", 30", 60" and 120" (i.e. 10, 20, 40 and 80 
h""'^ kpc at the distance of the cluster, respectively) are 
displayed in Figs. 7 to 10. 




Fig. 6. X-ray image obtained after removing the noise and the 
artefacts (see Appendix B). 



In the two first images (Figs. 7-8), we find significant 
(> 3(t) unresolved (i.e. not larger than one pixel) features. 
These generally have optical counterparts which would 
correspond to X-ray emitting galaxies. We find 12 galax- 
ies in our field (see Table 4). Our analysis is model inde- 
pendent, and allows the detection of X-ray sources. The 
method used by Dow & White (1995, a fitting procedure) 
gives actual measurements of the X-ray luminosities, but 
they only detect 4 galaxies in X-rays, and give upper lim- 
its for the others. Note that NGC 4874 and NGC 4889 are 
rather strong X-ray emitters (see Table 4). The strong X- 
ray source visible on Figs. 7-8 (labelled with a "Q") is a 
quasar. 

Structures are also observed around the positions of 
the two central galaxies at larger wavelet scales (see 
Figs. 9- 10). These detections clearly indicate the existence 
of two X-ray emitting groups around these two galaxies, in 
addition to the X-ray emission due to the galaxies them- 
selves. Vikhlinin et al. (1994) have already found these 
groups in the same ROSAT 'mia,ge, by applying a wavelet 
transform but at a single medium scale. 

The present analysis reveals two other features of par- 
ticular interest: 

1. A third maximum (labelled "A") which is clearly vis- 
ible in Figs. 8-9 between the two maxima associated 
with NGC 4874 and NGC 4889: this feature is not re- 
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Fig. 7. Wavelet image at a scale of 1 pixel (15"). The isophotes 
are drawn from 0.5 to 7 by steps of 0.5. The X and Y scales 
are drawn in units of ROSAT PSPC pixels of 15". The two 
brightest galaxies, NGC 4874 and NGC 4889, are indicated by 
x's. Other labels are: A for the secondary maximum, B for the 
filament between the two brightest galaxies, and Q for a quasar 
in the field. The positions of 10 bright galaxies are marked by 
crosses (see Table 2). 
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Fig. 8. Same as Fig. 7 at a scale of 2 pixels (30") 
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Fig. 9. Same as Fig. 7 at a scale of 4 pixels (60") 
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Fig. 10. Same as Fig. 7 at a scale of 8 pixels (120") 



lated to any particular galaxy and was also noted by 
White et al. (1993). It is located at GMP coordinates 
(-201,-54), which almost coincides with the position 
of the density maximum in the faint galaxy distribu- 
tion (the distance between the two peaks is only 5"). 
This X-ray maximum is indicated by an asterisk in the 
adaptive-kernel maps. 

A filamentary structure (labelled "B") between 
"A" and the galaxy NGC 4874 is visible in Figs. 8-9. 
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Table 4. Optical counterparts to X-ray features detected with 
the wavelet technique. Col. 1 : galaxy name; col. 2 : magnitude 
in the CMP catalogue; col. 3 : log Lx in the hard (0.4-2.4 keV) 
energy band from Dow & White (1995). 



Galaxy 


b 


log Lx 


name 






NGC 4874 


12.78 


42.19 


NGC 4889 


12.62 


42.11 


NGC 4898A 


14.85 


40.91 


NGC 4867 


15.44 


<40.91 


NGC 4883 


15.43 


<40.86 


IC 3959 


15.27 


40.79 


IC 3973 


15.32 


<40.79 


NGC 4869 


14.97 


<40.76 


NGC 4873 


15.15 


<40.69 


NGC 4864 


14.70 


<40.68 


NGC 4865 


14.54 


<40.62 


1257.3-1-2758 


15.15 


<40.47 



These two features can actually be seen directly on the 
reconstructed picture (Fig. 6). 

3.4. Comparison of optical and X-ray maps 

The distributions of the galaxies and X-ray gas are re- 
markably similar. In particular, the two brightest galax- 
ies are surrounded by galaxy groups that are also seen 
in the X-ray maps, and the faint galaxy population has 
a single density peak which is also found as a secondary 
maximum in the X-ray maps (feature "A"). The gas and 
galaxies seem therefore to be "swimming" in the same 
gravitational potential. However, there may exist features 
present on X-ray maps which are absent on optical maps. 
The response to the gravitational field of these two kinds of 
matter, which have very different intrinsic physical prop- 
erties, is different. Notice also that hydrodynamical sim- 
ulations predict emitting processes from the plasma not 
directly related to the potential well (Evrard 1990; Roet- 
tiger et al. 1993). As we discuss below (see § 6.1), feature 
"B" may be an example of such processes. 

4. The central region of Coma: kinematics 

As we have seen in § 3.2, the faint galaxies trace a 
much more regular single-peaked structure than the bright 
galaxies, which are mostly grouped around the two bright- 
est galaxies. We now make use of the velocity data for in- 
vestigating the kinematics of these two populations. For 
this purpose we use the Ciright,v and the Cjaint,v samples. 
As we have velocities only for 49 % of the galaxies in the 
C'faint,v sample, we have checked the results concerning 
the kinematics of the faint galaxy population on the sam- 
ple C faint, I9,v , which has 76 % of its galaxies with known 
redshifts. 



4.1. The velocity gradient 

The sample Ciright,v is dominated by the two galaxy 
groups around the central dominant galaxies; this can be 
seen by splitting the sample in two subsamples, one with 
velocities v < 6500 km s~"'^(and therefore closer to the ve- 
locity of NGC 4889), and the other with v > 6500 km s'^. 
There are 38 and 80 galaxies in the two subsamples respec- 
tively. The adaptive kernel density maps of these two sub- 
samples (which are 99 % magnitude complete) are shown 
in Fig. 11 (upper and lower panel, respectively, for the low 
and high velocity sample). 

On the other hand, the sample Cjaint,v does not show 
evidence for the presence of the two groups even if we con- 
sider separately the two subsamples at velocities respec- 
tively lower and higher than 6500 km s~^. Instead, there is 
evidence for a continuous velocity gradient in the opposite 
direction to the gradient between the two central galaxies; 
this is shown in the adaptive kernel map of the average 
velocity of this sample (Fig. 12; note that we have used 
an initial kernel size three times as large as that used in 
the other maps, i.e. 600" in this case, in order to increase 
the signal-to-noise in this map; see Appendix A for a more 
thorough explanation). Such a gradient has also been no- 
ticed by CD95 (see also Durret et al. 1995). Since this 
figure is drawn from a sample defined in a circular region, 
no preferential direction has been artificially introduced. 

Since the velocity field of the Cjaint,v sample of galax- 
ies is complex, we have looked for the direction along 
which the strongest correlation is obtained between the 
velocities and the coordinates along this direction, in or- 
der to have a quantitative estimate of the position an- 
gle and statistical significance of this gradient. The posi- 
tion angle of the maximal correlation is 40° (defined an- 
ticlockwise from North), with a rather large uncertainty 
of ±28°, estimated by performing 1000 bootstrap resam- 
plings of the data. This large uncertainty could be due 
to the existence of two superposed gradients, one in the 
East-West and the other in the North-South directions, 
as suggested by Fig. 12. The velocity gradient along the 
axis at position angle 40° is 0.25 ± 0.13 km s""'^ arcsec""'^ 
(375±195 km s""'^ h Mpc""'^), as deduced from a regression 
analysis between velocities and coordinates along this axis. 

The statistical significance of this gradient computed 
by the non-parametric Kendall correlation coefficient be- 
tween velocities and positions is 0.988. Its reality is further 
confirmed by the analysis of the more complete sample 
C faint, I9,v '■ the direction of the maximal correlation is al- 
most the same, 35° ± 23°, and the statistical significance 
is even higher, 0.998, despite the fact that the number of 
galaxies is lower in this sample; the gradient along the 
axis at position angle 35° is 0.47 ± 0.15 km s"-'^ arcsec"-'^ 
(705 ±225 km s~^ h Mpc~^). This gradient is much larger 
than the one found before, suggesting that incompleteness 
in the sample Cjaint,v has led to an underestimation of the 
true gradient intensity. 
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Fig. 11. Adaptive kernel density map of the galaxies in the 
sample Cbrtght,v with v < 6500 km s~^ (upper panel) and with 
V > 6500 km s~^ (lower panel). The contour spacing and the 
symbols are the same as in Fig. 1. 

Consistently with the lack of any coherent features in 
the average velocity maps of the sample and Ciright,v , 
we find no significant gradients on these samples, nor do 
we detect any significant gradient outside this inner circle 
of radius 1500", consistent with previous findings (Gre- 
gory 1975; Gregory & Tifl^t 1976). 

Fig. 13 shows the running velocity profile for galax- 
ies in the sample Cjaint,v vs. their positions along the di- 
rection of the velocity gradient (position angle 35°); the 
±1(7 confidence bands are estimated via 1000 bootstrap 
resamplings. The running velocity profile along the axis 




Fig. 12. Mean velocity adaptive kernel map for faint galaxies 
(C/asnt.tj sample). Contour levels are spaced by 100 km s""', 
and labelled with their mean velocity value. An initial kernel 
size three times as large as that used in the previous maps has 
been used here, in order to increase the signal-to-noise in this 
figure. Symbols are the same as in Fig. 1. 
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Fig. 13. Running velocity profile (heavy line) of the 
Cfaint,v sample along the gradient direction (35°); the dashed 
lines indicate the la upper and lower confidence bands as ob- 
tained via 1000 boostrap resamplings; the thin line in the mid- 
dle is the average of the boostrap resamplings. The two bright- 
est galaxies of Coma are indicated by squares; error-bars de- 
note the observational uncertainties on their velocities. 
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Table 5. Characteristics of the velocity distributions 

Sample mean dispersion skewness kurtosis Tail Index Omnibus-test 



km s ^ km s ^ probability 

TT, 6901 ± 86 1140 ± 56 -0.2 ± 0.2 0.4 ± 0.3 LOl (194 

CbTtght,v 6827 ±105 1119 ±81 0.0 ± 0.2 0.2 ± 0.4 0.93 0.59 

Cfatnt,v 6967 ± 102 1158 ±81 -0.4 ± 0.2 0.6 ± 0.4 0.98 0.98 

Caii,^ 6955 ± 56 972 ±35 -0.3 ± 0.1 0.4 ± 0.2 1.00 > 0.99 



at position angle of 40° as well as the running velocity 
profile of the C/aj„(_i9_t, sample are similar and are not 
displayed here. It must be noticed that neither NGC 4874 
nor NGC 4889 are located on the running velocity curve, 
and their relative velocity gradient has an opposite direc- 
tion to the gradient of the main cluster body. 

Velocity distributions 

The histograms and characteristics of the velocity samples 
C'v , C\right,v and Cfaint,v are displayed in Fig. 14 and Ta- 
ble 5. The values indicated in Table 5 are calculated with a 
biweight analysis (see, e.g.. Beers et al. 1990), and keeping 
into account the errors on velocities and the reduction to 
the cluster reference frame (Danese et al. 1980). All mean 
velocities and values of a are compatible with one another. 
The values we find for the mean cluster velocity and ve- 
locity dispersion are not different from those calculated by 
CD95. 

The velocity of NGC 4874 is 7176 ±15 km s'^, a value 
larger than the average cluster velocity by 275 km s""'^. 
The velocity of NGC 4889 is 6497 ± 19 km s'^, smaller 
than the mean velocity of Ct, by 404 km s""'^. These results 
are significant at levels larger than 4(7. The two central 
dominant galaxies are essentially in a symmetrical situ- 
ation (within K, 100 km s""'^) relatively to the Ct, mean 
velocity. 

In Table 5 we also list the values for the skewness, the 
kurtosis, and the Tail Index, relative to the Gaussian dis- 
tribution, and the probability of rejection of gaussianity 
according to an omnibus test, for the various samples. The 
Tail Index (see, e.g.. Beers et al. 1991) is a robust measure 
of the tails in a distribution, being lower than unity when 
the tails are underpopulated relative to a Gaussian. The 
omnibus test makes use of both the skewness and the kur- 
tosis to provide the probability that a given distribution 
is significantly not gaussian (see, e.g., D'Agostino 1986). 

Using a range of statistical tests, CD95 have shown 
that the velocity distribution of the whole GMP field, de- 
parts significantly from a gaussian. This is confirmed by 
our analysis (see Table 5, sample Caii^v )• On the other 
hand, the galaxy velocity distribution in the central clus- 
ter region (sample ) is only slightly non-gaussian. This 
non-gaussianity mainly arises from the faint galaxy popu- 
lation. In fact, the Ciright,v velocity distribution is not sig- 
nificantly different from a Gaussian, even if the low value 
of the Tail Index suggests the superposition of two popula- 
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Fig. 14. Histogram of velocities of the sample C„ (upper 
panel), Cbright,v (middle panel) and C faint, v (lower panel) be- 
fore (heavy-line histograms) and after (shaded histograms) the 
subtraction of the average velocity gradient from the observed 
velocities of galaxies in the Cfaint,v sample (which is also part 
of the C„ sample). The number of bins is approximately equal 
to the square-root of the number of galaxies in each sample. 



tions with different mean velocity in the data-set, probably 
the two groups around the central dominant galaxies (see 
Fig. 11 in § 4.1). On the other hand, the C/aj„(_t, velocity 
distribution is significantly non-gaussian, mainly because 
of its skewness. It is remarkable that if we subtract the 
average gradient (see § 4.1) from the observed velocities 
of this sample, the resulting velocity distribution is closer 
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Table 6. Properties of 13 BCMs and their "groups" 
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to a gaussian (the omnibus test rejects gaussianity with a 
probability of only 0.92). This is illustrated in Fig. 14. 

5. The outer regions of Coma 
5.1. Other groups 

Mellier et al. (1988) claimed the existence of a population 
of galaxies bound to several of the brightest members of 
the Coma cluster. We have repeated here their analysis 
using our larger data-set (see also Biviano et al. 1994). 
We have selected 13 galaxies with & < 15 (BCMs here- 
after, for Brightest Cluster Members) which are located 
near a density enhancement, as visible in an adaptive ker- 
nel map of the Call sample (Fig. 15). The names of these 
BCMs, positions, magnitudes and velocities are listed in 
Table 6. While some of the peaks visible in Fig. 15 may 
not correspond to real structures, the purpose here is to 
select all possible groups, and to run a statistical test for 
the presence of two populations within these groups. 

We have constructed 13 "groups" centered on these 
BCMs by taking all the galaxies closer than 250" from each 
BCM; the radius of 250" has been chosen as a compro- 
mise between the need of selecting galaxies close enough 
to the BCMs, and still having sufficient statistics. Note 
that with such a radius, the "groups" of NGC 4874 and 
NGC 4889 would be superposed; we have chosen to assign 
to the NGC 4874 "group" the galaxies falling in the su- 
perposition zone, since from Fig. 4, NGC 4874 seems to sit 
on a larger density peak. 

In Table 6 we list in col. (7) the ratio between the num- 
ber of galaxies with available velocities, N^, and the total 
number of galaxies, N, in each "group", BCMs excluded. 
Since we do not have any galaxy with available velocity 
in the vicinity of 3 BCMs, our analysis is effectively re- 
stricted to 10 putative "groups". We also give the aver- 
age velocities for 4 "groups" with enough galaxy velocities 
available; these average velocities (relative to the velocity 
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Fig. 15. Adaptive kernel density map of all the cluster mem- 
bers in the GMP region (sample Call )• The x's indicate all 
cluster members brighter than 6 = 15 in this field. Contour 
spacing is the same as in Fig. 1, but contours in the central 
region have been omitted for clarity. 



of the BCM) are listed in col. (8) of Table 6 (these aver- 
ages are computed without including BCMs). The average 
"group" velocities are equal to the BCM velocities, within 
the errors; however, these errors are quite large. In order to 
check for the existence of a population of galaxies bound 
to the BCMs, we have merged together the velocities (rel- 
ative to the velocities of the BCMs) of the 74 galaxies 
assigned to the 10 "groups" (BCMs excluded). The null 
hypothesis that we are testing is that this distribution is 
the same as that of all the cluster galaxy velocities relative 
to the cluster mean. This should be the case, if the sample 
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of galaxies in the vicinity of BCMs is a random sample of 
the cluster population. 

A Kolmogorov-Smirnov test (see e.g. Stephens 1986) 
does not reject the hypothesis that the 74 "group" galaxy 
velocities (relative to their BCM velocities) are drawn 
from the same parent distribution of the 469 velocities (rel- 
ative to the cluster average velocity) in the sample Caii^v ■ 
However, there is a difference between the velocity distri- 
butions of the "group" galaxies and of the Cau,v galaxies, 
when only galaxies with h < 15.5 are considered (there 
are 13 and 89 galaxies brighter than this magnitude, in 
the "group" and the Cau,v sample respectively). This dif- 
ference is significant at the 94.01 % confidence level, ac- 
cording to the Kolmogorov-Smirnov test. In fact, the ve- 
locity dispersion of the 13 "group" galaxies brighter than 
h = 15.5, with respect to the velocities of their BCMs, is 
only 424 km s~^. We conclude that the population bound 
to the BCMs, discovered by Mellier et al. (1988), is made 
of bright galaxies only. A similar result was found in § 3.2 
for the central region of Coma: bright galaxies tend to lie 
in sub clusters and faint galaxies trace the overall cluster 
potential. 



5.2. The region of NGC 4839 

Let us now concentrate on the well known group around 
the giant galaxy NGC 4839. The presence of this group 
makes the velocity distribution of the Cau,v sample signif- 
icantly skewed and therefore non-gaussian (see Table 5). 
This was recently pointed out by CD95, who were able 
to split the group and the cluster components using the 
KMM algorithm (see, e.g.. Bird 1994); they obtained an 
average velocity of 7339 km s~^ and a velocity dispersion 
of 329 km s~^ for the group. The average velocity is in per- 
fect agreement with the velocity of NGC 4839 (7317±20 
km s""'^). 

However, CD95 did not address the possible problem 
of the incompleteness of their sample; in order to check 
the robustness of their result, we have considered only the 
galaxies of sample Cau,v with magnitudes h < 17; the com- 
pleteness of this sample is 82 % (see Table 1). By perform- 
ing a maximum likelihood fit with two gaussians to the ve- 
locity distribution of this sample, we recover (within the 
errors) the average velocity and velocity dispersion found 
by CD95 for the group population. The nice agreement 
with the values found by CD95 suggests that their result 
is robust, despite incompleteness. Using the ratio between 
the numbers of galaxies belonging to the two gaussians 
and assuming that all the likely cluster members with- 
out available velocities belong to the cluster and not to 
the group, we can set a lower limit of 0.17 to the ratio 
between the number of galaxies in the group and in the 
cluster. 



5.3. The distribution of star-forming galaxies 

As noted by many authors (Sodre et al. 1989; Zabludoff 
& Franx 1993; CD95; Donas et al. 1995; Gavazzi et al. 
1995), the star-forming and the quiescent galaxy popula- 
tions in Coma have markedly different velocity distribu- 
tions. We have separated these two populations by their 
colours, respectively larger or equal to, and smaller than 
h — r = 1.7. In the sample Ca/; there are 278 blue clus- 
ter members, of which 64 with measured velocities, and 
1073 red cluster members, of which 384 with measured 
velocities. The density profiles and velocity distributions 
of the red and blue Coma cluster members are shown in 
Fig. 16. The density profile of the blue galaxies is fiat- 
ter than that of the red galaxies in the center, yet the 
two profiles are similar beyond 1200" (corresponding to 
0.8 h~^ Mpc) from the center; the velocity distributions 
of the two populations are different (0.991 probability, ac- 
cording to a Kolmogorov-Smirnov test), both because of 
a different mean and a different dispersion (the velocity 
dispersion of the blue galaxies is 30 % larger than that of 
the red galaxies). 

The sample Cau,v being rather incomplete, we have 
performed a similar analysis on the 82 % complete sample 
of galaxies brighter than h < 17, and found approximately 
the same results. 

Caldwell et al. (1993) have identified 17 early-type 
galaxies with an abnormal spectrum in the Coma clus- 
ter. Given the peculiarity of their spectrum, most of these 
galaxies are thought to have suffered a starburst episode 
~ 1 Gyr ago (Caldwell et al. 1995). A possible explanation 
for such a starburst is a recent interaction with the clus- 
ter environment. These 17 galaxies share a similar veloc- 
ity distribution as that of the blue galaxies (see Fig. 16). 
Their average velocity (7526±92 km s""'^) and velocity dis- 
persion (297±1050 km s"-'^) are also quite similar to those 
of the NGC 4839 group, although their spatial distribution 
is not restricted to the region of NGC 4839. Note that the 
values of the average velocity and dispersion are different 
from those given in Caldwell et al. (1993) and in CD95, 
because here we have used the biweight estimator, which 
is more reliable than the classical dispersion in the case of 
poor data-sets (Beers et al. 1990). While the average ve- 
locity is roughly two standard deviations offset from the 
average velocity of the NGC 4839 group, and the velocity 
dispersion is affected by an enormous error, this result is 
intriguing. 

6. Discussion 

6.L The mam cluster body 

As shown in the previous sections, the bright and faint 
galaxies have markedly different distributions and kine- 
matics. Bright galaxies tend to be located in groups domi- 
nated by their very brightest members, while faint galaxies 
delineate a much smoother single-peaked structure. This is 
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Fig. 16. Upper panel: the adaptive kernel number density pro- 
files of blue (thick line) and red (thin line) cluster members. 
The dashed lines indicate the 1-a confidence bands on the 
two profiles, obtained with 1000 bootstrap resamplings. Lower 
panel: the velocity distribution of the blue (thick histogram) 
and red (thin histogram) galaxy populations, and of the 17 
galaxies with abnormal spectra (dashed histogram) detected 
by Caldwell et al. (1993). 

true not only for the groups around the two core dominant 
galaxies, but also for other groups even in the external re- 
gions. This is reminiscent of the luminosity segregation 
already observed by Rood (1969). The dichotomy is nat- 
urally not perfect, and there probably are faint galaxies 
located in groups as well, as suggested by the Western ex- 
tension (towards NGC 4874) of the maximum in the faint 
galaxy density map (Fig. 5). The X-ray maps show sev- 
eral maxima; in particular, two are centered on NGC 4889 
and 4874, and one is coincident with the faint galaxy dis- 
tribution density peak. The coincidence of the optical and 
X-ray emitting gas distributions strongly suggests that the 
detected features are real. 

It is usually assumed that the main cluster structure 
coincides with the overdensity around NGC 4874. How- 
ever, neither this galaxy nor the brightest cluster member 
(NGC 4889) are at rest in the cluster potential, since their 



velocities are different from the average cluster velocity; 
they are more likely to be located in subclusters still in 
the process of merging with the cluster. Indeed, the den- 
sity and X-ray maxima around the two central dominant 
galaxies indicate the presence of groups linked to these 
two galaxies; this is supported by the kinematical analysis 
(see Fig. 11), and by previous investigations (Fitchett & 
Webster 1987; White et al. 1993; Vikhlinin et al. 1994). 
CD95 claimed that the two central dominant galaxies do 
not lie at the mean velocity of their own groups; but they 
did not realize that the bright galaxies do cluster around 
the two central galaxies, while only the faint galaxies do 
not (indeed their velocity gradient is opposite to the mean 
velocity gradient of these two groups, see § 4.1 and § 6.2). 

The X-ray secondary maximum (feature "A" in Fig. 9) 
has already been noted by White et al. (1993), who pro- 
posed to identify it with a remnant of gas stripped from 
NGC 4889. The fact that this peak is also seen in the 
(faint) galaxy distribution suggests an alternative expla- 
nation: this could be the X-ray maximum of the mam body 
of the Coma cluster. The "A" feature indeed possesses 
another property typical of clusters: it is visible on small 
or medium scales, while it is blurred on the largest scale, 
indicating that its intrinsic characteristic scale is small, 
i.e. ~ 60" (corresponding to 40 h~^ kpc). Indeed, in two 
papers devoted to the analysis of 14 X-ray clusters, with 
improved techniques, Gerbal et al. (1992) and Durret et 
al. (1994) have shown that X-ray core radii are very small, 
~ 50-100 h~^ kpc (in opposition to what is often claimed 
in the literature, but in agreement with results derived 
from the gravitational lensing analysis of the unseen mat- 
ter distribution). The central peak of the X-ray brightness 
must then have a medium typical scale, and this is why 
it is visible on the medium scale image, and not on the 
larger one. 

Our interpretation naturally accounts for the fact that 
both NGC 4889, and NGC 4874 have velocities signifi- 
cantly different from the cluster mean (see § 4.2), since 
both galaxies are located within their own groups, even- 
tually falling on to the main body of Coma. These groups 
would then only be two of the groups merging with, or 
falling onto the main cluster body. In this case. Coma 
would still be undergoing a formation process, with the 
main body continuously accreting matter, i.e. more or less 
large groups of galaxies. The main body would in fact be 
formed by the mixing of such groups and clusters arrived 
in the past, suggesting the continuous existence of violent 
relaxation processes. Notice that such processes are noth- 
ing but secondary mfall, adapted to a more realistic view 
of the universe than the simple case of spherical symmetry. 

In particular, we suggest that the infall of the 
NGC 4874 group onto the high density cluster core (region 
"A" in Fig. 9), is producing a compression wave that has 
increased the pressure, and, as a consequence, the X-ray 
luminosity of the gas along the east edge of this group; the 
X-ray emission of this compressed gas would correspond 
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to the extended filament "B" in Fig. 9. Such an interpre- 
tation finds support in hydrodynamical simulations (see 
e.g. Evrard, 1990) showing that the X-ray emission is not 
only associated with a smooth gas distribution, but can 
also arise from "shocks" induced by the arrival of groups 
into the cluster. 

Finally, the fact that we detect the main body in 
the faint galaxy population, while the bright galaxies are 
mainly located in subclusters, can be explained as the ef- 
fect of dynamical evolution in the groups before their infall 
onto the cluster. The very existence of such large galaxies 
like NGC 4889 and NGC 4874 argues in favour of this hy- 
pothesis. In fact, brightest cluster members are expected 
to form in groups before violent relaxation occurs (Mer- 
rill 1984), probably through mergers which imply small 
velocity encounters. Such small velocities may be attained 
among bright, massive galaxies by the loss of their kinetic 
energy through the process of dynamical friction (see, e.g., 
Biviano et al. 1992; Yepes & Dominguez-Tenreiro 1992). 
The groups could then evolve to a core-halo structure, 
with bright galaxies in their core, and faint galaxies in 
their halo. When the infall onto the cluster occurs, the 
less bound halo population would be more easily tidally 
stripped off the group, and become part of the cluster 
population, at variance with the core population which 
could resist longer to tidal disruption. Of course there are 
probably unevolved groups as well infalling onto the clus- 
ter, but these would loose their identity very soon because 
of tidal disruption, and escape detection; the dense cores 
of evolved groups would survive longer, as indicated by 
numerical simulations (Gonzalez-Casado et al. 1994). 

6.2. Does the mam body rotate? 

As shown in § 4.1, the sample C/aj„(_t, has kinematical 
properties on its own; this lends further support to the 
hypothesis that the density peak of this sample is inde- 
pendent from the density peaks associated with the two 
central dominant galaxies. In particular, we have shown 
that the faint galaxy velocity distribution has a significant 
gradient (see Fig. 12 and 13). This gradient was also noted 
by CD95, but they did not recognize that it is restricted 
to the faint galaxy population (no significant gradient is 
found in the velocity distribution of the bright galaxy pop- 
ulation; if anything the gradient is in the opposite sense, 
see Fig. 11). 

It is quite remarkable that this velocity gradient is 
globally oriented from the North-East to the South- West 
(35°-40°), in rough coincidence with several main direc- 
tions in Coma. These privileged directions have recently 
been pointed out by West (1994), who explains them as 
a consequence of the infall of structures onto the cluster 
from the large scale structure filaments surrounding it. In 
particular, the direction of the velocity gradient is closer to 
the direction between the centre of Coma and NGC 4839, 
rather than to the direction to the neighboring cluster 



Abell 1367. Note, however, that the NGC 4839 group av- 
erage velocity (~ 7400 km s""'^) is opposite to what ex- 
pected from an extrapolation of the central gradient to 
the group position. 

This gradient is an intrinsic property of the main clus- 
ter body. In fact, it cannot simply refiect the presence of 
sub clustering around the two central dominant galaxies, 
because it goes in the opposite sense. On the other hand, 
the location of the Coma cluster at the intersection of 
several sheets may suggest that the velocity gradient re- 
fiects the large-scale-structure galaxy distribution. In par- 
ticular. Coma is within the Great Wall (Geller & Huchra 
1989; Ramella et al. 1992), which is inclined with respect 
to the line of sight. However, while in the same direction, 
the intensity of the gradient found here is much higher 
than that of the Great Wall, which is indeed roughly sim- 
ilar to the one previously detected by Gregory & Tifft 
(1976; ~ 0.1 km s"-'^ arcsec"-'^ , i.e. ~ 15 km s"-'^ h Mpc"-'^). 
Morevoer, we only detect a gradient in the inner circle of 
1500" radius. 

What is the physical meaning of this gradient? The 
velocity map (see Fig. 12) is complex and not as ex- 
pected from simple rotation. Moreover, if we interpret the 
gradient as arising from rotation, the rotational velocity 
would be too low to account for the fiattening of the main 
body, as one can estimate using the anisotropy parameter 
(v/a)* defined in, e.g., Ferguson & Binggeli (1994, eq.5). 
The fact that the gradient is principally oriented towards 
NGC 4839 may suggest a transient tidal effect. Among 
the numerous consequences of such a transient tidal ef- 
fect, are the deformation of the matter distribution, and 
a mean velocity field which is no longer constant. The ef- 
fects of tides on spheroidal galaxies infalling onto our own 
Galaxy are described in a paper by Piatek & Pryor (1995). 
They find that tides deform the spheroidal galaxies and 
induce apparent rotation: tides produce large ordered mo- 
tions rather than large random motions. We think that 
something similar may be happening in the Coma cluster, 
with the central part of the main body suffering from the 
tidal infiuence of the infalling groups, like those around 
NGC 4874 and NGC 4839. The search for this kind of ef- 
fect with realistic numerical simulations of clusters could 
give constraints on the relative masses of the satellite and 
main body, as well as on the history of encounters. 

6.3. Ongoing mfall 

The groups around NGC 4874 and NGC 4889 are likely 
to be orbiting within the Coma cluster, and they are not 
at rest at the bottom of the cluster potential. There are 
at least two more groups, one centered on the galaxy 
NGC 4839 and the other, less conspicuous, on the galaxy 
NGC 4911, and other groups maybe present around other 
bright galaxies (Mellier et al. 1988). The average velocities 
of these galaxies and their groups are also quite different 
from the cluster mean velocity. Altogether, this suggests 
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that the Coma cluster is still in the process of accreting 
groups from the surrounding large scale structure. 

Not only is Coma accreting groups but also individual 
galaxies. The flatness of the density profile of the blue 
galaxy population near the center indicates that these 
galaxies tend to avoid the center, yet the similarity of their 
density profile to that of red galaxies (see Fig. 16) at large 
radii suggests that these blue galaxies are cluster mem- 
bers. It is then possible that these galaxies spend only a 
little time in the central region of the cluster, because they 
have very elongated orbits. This is also indicated by the 
double peaked velocity histogram of the blue galaxy pop- 
ulation, typical of galaxies on strongly radial orbits (see, 
e.g., Merritt 1987). This situation is not at all unusual in 
galaxy clusters (Moss & Dickens 1977; Sodre et al. 1989; 
Zabludoff & Franx 1993; Biviano et al. 1995c). Supporting 
this interpretation is the similarity of the blue galaxy ve- 
locity distribution to those of the Hl-deficient galaxies (see 
Fig. 3 in Gavazzi 1987), and of the abnormal-spectrum 
early-type galaxies (see Fig. 16) detected by Caldwell et al. 
(1993). Both the Hl-deficient galaxies and the abnormal- 
spectrum early-type galaxies are likely to have interacted 
with the cluster environment, and are therefore probable 
cluster members, so that this may also be true of the gen- 
eral blue galaxy population. 

The group around NGC 4839 has an average veloc- 
ity close to the peak in the velocity distribution of the 
abnormal-spectrum galaxies. This fact led CD95 to con- 
clude that both the group and these galaxies are falling 
onto the cluster coming from the same galaxy filament. 
While this is certainly a possibility, supported by the 
similarity of the velocity distributions of the abnormal- 
spectrum and blue galaxies, the abnormal-spectrum galax- 
ies also have a similar velocity dispersion as the NGC 4839 
group, and most of them are located near this group. The 
possibility then remains that these galaxies are (or have 
been) members of the NGC 4839 group. In this case, if 
they have already interacted with the cluster, this would 
be true of the NGC 4839 group as well. 

Burns et al. (1994) suggested that the NGC 4839 group 
has passed through the core of Coma following a straight 
trajectory. As shown by Gonzalez-Casado et al. (1994), 
this would cause the group structure to be severely af- 
fected by tidal forces. However, the NGC 4839 group has 
a richness, velocity dispersion, and X-ray luminosity that 
fits very well the relations between these quantities, as 
derived for clusters of galaxies (Edge & Stewart 1991; Gi- 
rardi et al. 1993); this argues against the possibility that 
the group has suffered strong tidal effects. A possible al- 
ternative scenario is a non-radial infall of the NGC 4839 
group onto the cluster; one or several passages through 
the outer parts of the cluster may already have occurred, 
which would leave the group intact, except for some galax- 
ies that may have been lost in its trail. Such a scenario 
agrees with predictions derived from recent numerical sim- 



ulations in which the effects of dynamical friction are in- 
cluded (see e.g. Chan et al. 1995 and references therein). 

7. Conclusions 

The analysis of Coma simultaneously in X-rays and at 
optical wavelengths allows us to draw a coherent picture 
of this cluster: 

1. A main body centered on the peak of the distribution 
of faint galaxies, which coincides with a secondary X- 
ray peak, is unveiled; groups are falling on to this main 
body; 

2. Faint galaxies are better tracers of the cluster struc- 
ture than bright galaxies, since many of the latter are 
subclustered around the very brightest galaxies; 

3. We detect a gradient in the velocities of the faint galax- 
ies tracing the main cluster body; we suggest that 
this gradient arises from tidal interactions of the main 
cluster body with infalling groups, in particular with 
NGC 4839; 

4. The infall of the group surrounding NGC 4874 may 
have induced a pressure wave in the X-ray emitting 
gas; 

5. Groups are still visible mainly because of a bound pop- 
ulation of bright galaxies, which has probably evolved 
to a dense core structure before the infall of their own 
group, and is thus less prone to tidal stripping; 

6. The cluster is still accreting these groups and also indi- 
vidual galaxies from the surrounding large scale struc- 
ture. 

Acknowledgements. We deeply thank J.G. Godwin, N. Met- 
calfe & J.V. Peach for providing us with their unpublished 
catalogue of stellar objects in the Coma field, M. CoUess for 
providing us with his catalogue of redshifts in advance of publi- 
cation, and C. Clemens for sending us an electronic copy of the 
ZCAT catalogue. We are also grateful to S. White for providing 
us with the /SO^AT images of Coma before their release in the 
Garching public archive. We acknowledge the help of D. Fadda 
in the software implementation of the adaptive kernel tech- 
nique at an early stage of this project. We thank D. Merritt 
and A. Pisani for useful discussions. We acknowledge financial 
support from GDR "Cosmologie" , CNRS, and from the Euro- 
pean Cosmology Network (a E.U. H.C.M program). C. Lobo 
is fully supported by the BD/2772/93RM grant attributed by 
JNICT, Portugal. 

Appendix A: Adaptive kernel maps 

A thorough description of the adaptive kernel techniques is 
given in Pisani (1993, 1995). Here we only address the issues 
that are not reported in those papers. 

A.l. Statistical significance 

In order to assign the statistical significance to a given fea- 
ture we have run 1000 bootstrap resamplings on each map. In 
practice, this consists in evaluating an adaptive kernel map of 
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the quantity of interest (typically, the density) for each boot- 
strap sample; we then compute the standard deviation of the 
bootstrap values for each pixel in the map. This provides an 
estimate of the error in the observed quantity in each pixel. 
We then subtract to the observed value in each pixel a given 
number (say, 3) of standard deviations, and obtain a new map 
which we compare to the original one, to see if the observed 
features are still visible. Features affected by large errors will 
not be visible in the density map at —3(7 confidence level. 

A. 2. Velocity maps 

We define an adaptive kernel map of the average velocities 
by multiplying the 2-dimensional gaussian adaptive kernel, 
K^^^\ x ), of each galaxy j, by its velocity (when available), 

v{i)=Y.Kr'(^>^iT.^f"' (Al) 

where x is the position vector of a given pixel of the image, 
and is the number of galaxies with available velocity. 

In order to increase the signal-to-noise in these maps, we 
adopt a larger initial kernel size than in the density maps. 

A. 3. Substructure maps 

Before defining a 5-map, let us remind the definition of the 
parameter 8: 

^3 = + (''■f.j ~ '^^ f (A2) 

where Vj is the average velocity and a^^j is the velocity disper- 
sion of the galaxy group defined by the a/TVi, nearest neigh- 
bours to the j-th galaxy, and v and are the average velocity 
and velocity dispersion of all the galaxies with available 
velocities. 

We then define a local average velocity v(x), as in eq. (Al), 
and a local velocity dispersion, (7„(x), as follows: 

al(x) = Y^Kf^\x)(v,-v(x)flY,Kf'''^ (A3) 

Then, as a straightforward extension of the Dressier & 
Shectman (1988) test, we compute the local values of 8: 

S(i) = {v(i)-vf Jr{<yv(i)-<yvf (A4) 
The isopleths of 8(x) give the 5-map. 

Appendix B: Wavelet maps 

Images which exhibit features of different characteristic sizes 
at various locations are much more accurately analyzed using 
space-scale techniques than monoscale informations. In Slezak 
et al. (1994), we discussed the ability of the wavelet transform 
to locate these structures according to their typical scales and 
presented quite extensively the basic mathematical formulae 
as well as a fast computation algorithm. Therefore, we focus 
especially this appendix on the use of a multiscale vision model 



for restoring an image where some components have been sub- 
tracted. We have already applied these concepts to design a 
hierarchical filter which allows to remove the small-scale noise 
while keeping the relevant details associated to the objects 
(Slezak et al. 1994). Both improvements and generalization of 
this method are explained here in a concise way since a full 
description of the model we use and numerical simulations and 
tests will be given in Rue & Bijaoui (1995). 

B.l. The multiscale vision model 

Using an isotropic analyzing wavelet without any decimation 
from one scale a, = 2' to another, a 2D image is transformed 
into a set of wavelet images {VUj} exhibiting details related to 
different spatial frequencies, so that a complex object showing 
features of different sizes is associated to non zero coefficients 
at successive scales. An object now has to be described in this 
wavelet space and by consequence the connectivity property 
which defines it in the direct space needs to be enlarged to the 
3D set of wavelet coefficients W,{k, I). To do so, one must define 
an interscale connectivity relation linking the various compo- 
nents corresponding to the connected regions D,{n) detected 
at each scale a, by the field labelling L,{k,l) = n. This image 
segmentation permits to select significant wavelet coefficients 
from noisy data. Rue & Bijaoui (1995) decided to base this 
relationship TZ on the location (kmax, Imax) of the maximum 
value inside D,{n), and we have adopted the same rule : 

D,{n) Tl D,^l{m) if ij + l (fcmax, ^max) = m . (Bl) 

Starting from the largest scale to the smallest one, the iden- 
tification of the fields which appear to be connected in this 
manner results in a connection tree. So, any complex object 
is viewed as an arborescent structure V in the wavelet space, 
the whole information of which is required to characterize cor- 
rectly the object with all its components. A key issue in this 
approach is thus the possibility to extract easily these compo- 
nents by considering their related sub-trees, i.e. the regions of 
the wavelet structure connected to the field identified as the 
largest structural element of the component. One can there- 
fore either study these well-defined components in a separate 
way or analyze an object where undesirable features like su- 
perimposed smaller objects have been quite perfectly removed 
by eliminating their sub-trees. Note that the partition related 
to such somewhat bright sub-objects can be distinguished from 
genuine components belonging to a complex object by the max- 
imum wavelet value inside their root domain. It must be higher 
than the values of the maxima for the connected fields at the 
previous larger and next smaller scales since wavelet values in- 
deed increase or decrease according to the match between the 
wavelet scale and the typical size and amplitude of the signal. 
The most efficient way to characterize each source detected in 
the wavelet space is to analyze an image / of this source in the 
direct space. Let us now address this issue. 

B.2. The image reconstruction 

Assuming a Gaussian noise and providing that the standard de- 
viation parameters a, of the Gaussian density probability dis- 
tributions for the noise process, P{W,{k, I)), are obtained from 
the clipped histogram of the oversampled and noise-dominated 
coefficient map Wi (see Slezak et al. 1994), the fields D,{n) are 
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labelled by performing an image segmentation at a 3 (Tj level 
on each wavelet map . The connection trees are then con- 
structed following the rule described in § B.l. They are pruned 
in order to remove isolated domains and the existing objects 
and components are identified from the maximum value of their 
root fields. 

We now have to compute for each object an image / from the 
set V of the wavelet coefficients associated to its tree by solving 
the equation : 

V = O (/) = (Pt, WT) (/) = Wt> , (B2) 

where WT and Pd design the operators associated to the 
wavelet transform and to the projection in the domains "D = 
y^ ^. respectively. The mere addition of the thresholded 
frames is in fact unable to restore an image which will give 
again the same wavelet coefficients since the thresholding re- 
sults in a wavelet structure which cannot be considered as the 
wavelet transform of an image. 

This inverse problem has at least one solution: that corre- 
sponding to the image which has been processed. Among the 
class of mathematically acceptable solutions, the choice of the 
physically correct one is controlled by a regularization condi- 
tion. Many are available. For instance we first decided to take 
into account the consistency between V and W-p, and we used 
a simple direct algorithm similar to the iterative approach de- 
scribed in Van Cittert (1931) for which the regularization is 
achieved by constraining the support of the solution / (Slezak 
et al. 1994). Another choice is to minimize the energy of the 
residuals between the solution and the initial set of wavelet 
coefficients, i.e. 

\\V-0 (/)|| is minimum with ||V|| = Ey W^{k,l) , (B3) 

which leads to the so-called conjugate gradient algorithm. This 
method appears better suited for accurate image restoration 
of partial trees and of sub-trees than our previous technique 
designed mainly for restoring whole trees, and we decided to 
apply it to the ROSAT X-ray images of the Coma cluster. It 
is fully described in Rue & Bijaoui (1995), and only the key 
formulae of their iterative algorithm are given hereafter. 

Let us note O the joint operator which transforms a wavelet 
structure into an image and and FF^"'' the wavelet transform 
inside X* of / at step n. The joint operator O appears to be ba- 
sically equal to the addition of the wavelet frames [FFT(FF8)]8. 
We get : 

/("+!) = /(") _^ „(")/(") (B4) 

with = 6(V-PF4"^) + , where f"'' stands for 

the residual image at step n and a'-"-' and /?'•"-' are convergence 
parameters. At each step, the solution /("+^) is thresholded 
in order to get a positive image and the energy of the resid- 
ual wavelet coefficients V — 0(1^"'^^^) is tested against a given 
precision threshold. Initializations are performed by using the 
classical wavelet reconstruction operator adding all the wavelet 
planes to the smoothest approximation of the signal for 1^°^ 
and by setting /?'•''-' to zero. 

The reconstruction appears optimal only if the maximum 
of the wavelet coefficients involved for a given object is included 
into V. If it is not the case, an important information about 



the object is missing and the discrepancy between its original 
and restored images may be large. 

Thanks to the unfolding of the data in the wavelet space 
and by considering whole hierarchical connection trees or part 
of them, this iterative process enables one to restore either 
images of objects where the noise and superimposed features 
have been accurately subtracted without loss of any small-scale 
information, or images of these features where the extended 
underlying main objects have been removed. 
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